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Abstract 

Effects associated with the existence of isolated zero modes of Major ana fermions 
are discussed. It is argued that the quantization of this system necessarily contains 
highly extended quantum states and that populating and depopulating such states by 
interacting with the quantum system leads to long-ranged teleportation-like processes. 
Also leads to spontaneous violation of fermion parity symmetry. A quasi-realistic 
model consisting of a quantum wire embedded in a p-wave superconductor is discussed 
as an explicit example of a physical system with an isolated Majorana zero mode. 



1 Introduction 

It is a great pleasure to dedicate this article to the 100'th anniversary of the birth of 
Ettore Majorana. As a testimony to his lasting influence on science, we shall describe 
how one of his great insights, used in a modern context, can be related to a particular 
macroscopic quantum phenomenon. 

The idea is related to the observation by Majorana that a relativistic fermion 
such as the electron can be meaningfully decomposed into more basic degrees of free- 
dom, essentially by taking the real and imaginary parts of its wave- function pQ. In 
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relativistic field theory, what one obtains are called Majorana fermions, which have 
become the basic building blocks of super symmetric field theories and supply a sce- 
nario whereby the neutrinos which are observed in nature can have mass. We shall 
use this idea in a different context. In quantum condensed matter, the appearance of 
an emergent Majorana fermion would provide an excitation of a system that has min- 
imal degrees of freedom. The wave-function of the single-particle state would obey 
a Majorana condition, which would forbid quantum fluctuations of its phase. The 
utility of this fact has already been recognized in the context of quantum comput- 
ing |2j - In the present manuscript, we will elaborate on our previous observation 
[5] that in some cases this can provide isolated states with wave-functions which are 
peaked at multiple, well separated locations. In a controlled setting, this can be used 
to create a condensed matter realization of the Einstein-Podolsky-Rosen effect and 
even a version of teleportation by long-ranged tunnelling. 

Major ana's original motivation for inventing the Majorana fermion was to avoid 
the negative energy states that relativistic particles invariably seem to possess by 
identifying the negative and positive energy states of a relativistic wave equation as 
manifestations of the same quantum excitation. 

In second quantization, the positive energy state can be occupied by a particle. 
Filling a positive energy state creates an excited state of the system with positive 
energy. On the other hand, a negative energy state should be regarded as typically 
being already filled by a particle. An excitation of the system is then found by 
emptying the negative energy state, or creating a hole. The system is put in a higher 
energy state by removing a negative energy particle, equivalently, creating a positive 
energy hole. 

Majorana's idea can be implemented when there is a particle-hole symmetry. 
Then, for a given particle state, there exists a hole state with the same energy and 
with a wave-function that is related to the particle wave-function by a simple trans- 
formation. Then, by making the appropriate identification, one could indeed identify 
these as one and the same quantum state. Of course, the resulting system has half 
as many degrees of freedom. 1 

To illustrate the idea, let us recall the conventional second quantization of com- 
plex fermions, which could be either relativistic or non-relativistic. We begin with 
the assumption that in some approximation it makes sense to discuss a single non- 
interacting particle whose wave-function obeys the Schrodinger equation 

ih^{x,t) =H V(x,t) (1) 

where Hq is the single-particle Hamiltonian operator. Generally, as in the case of 
the Dirac equation, the Hamiltonian Hq could be a matrix, as well as a differential 
operator, and *ff(x, t) a column vector whose indices we shall suppress. 2 The second- 

1 For a comprehensive account of issues to do with positive and negative energy modes of relativistic 
bosons and fermions, see the series of papers |h]-|llj. 

2 An example is the Dirac Hamiltonian in 3+1-dimensions 

Hq = id ■ V + (3m 



where a and (3 are a set of four Hermitian, anti-commuting 4x4 Dirac matrices. There exists a matrix T 



quantized field operator also typically obeys this wave equation plus the equal-time 
anti-commutation relation 

{*(x,t),&(y,t)}=6(x-g) (2) 

It is this anti-commutator which defines *$>(x,t) as an operator. It can further be 
used to derive the wave equation (^Q) from the second quantized Hamiltonian, 

using the Hamilton equation of motion 

ih^*(x,t) = [V(x,t),H] 

We shall assume that Hq is a Hermitian operator which has eigenfunctions and a 
spectrum of real eigenvalues 

H iPe(x) = E%I)e{x) 

The energy E can be both positive and negative, in fact for the relativistic elec- 
tron, if Q were the Dirac equation, there are necessarily negative eigenvalues and 
the spectrum is unbounded below. The eigenfunctions obey the orthogonality and 
completeness relations 

/ dxil) ] E (x)^ E ,(x) = 8 E E' , ^2ipE(x)ip E (y) = 5(x - y) (4) 

J E 

The delta function and summation in these formulae should be understood in a 
generalized sense where they are a Kronecker delta and a sum for discrete components 
of the spectrum and a Dirac delta function and integral for continuum spectrum. 

In this system, one then forms the second quantized field operator by superposing 
the wave-functions with creation and annihilation operators, 

tf(x,t) = 4>E(x)e-' lEt / h a E + ^E(x)e- lEt ' h bl E 

E>0 E<0 

Here, a E is the annihilation operator for a particle with energy E and b_ E is the 
creation operator for a hole with energy —E. When they obey the algebra with 
non- vanishing anti-commutators 

{a E ,a E ,^ = 8 EE > , h- E ,tf_ E ,\ = 8 EE ! 



with the property Ta T = a* and T8T = — j3* , so that, THqT = —Hq and Tip^ = ip-E- This is a one-to-one 
mapping of positive to negative energy states. Explicitly, if the matrices are represented by a. = ^ 

with a the Pauli matrices and B = ( ? . | , then we can form the matrix T = ( , l* 7 I • Note that, 

yi OJ \ia* J 

in this case T = V* and T 2 = 1. A Majorana fermion obeys the reality condition ^(x,t) = T^!*(x, t). 



the $?(x, t) obeys the anticommutator (J2J). The completeness condition in Eq. (J5J) is 
essential for establishing this. 

The ground state of the system, |0 >, is the state where all positive energy levels 
are empty and where all negative energy levels are filled, or alternatively all hole states 
are empty. In the second quantized language, it is annihilated by the annihilation 
operators, 

a E \0 >= = b_ E \0 > 

Excited states are created by operating on |0 > with a E and b_ E . The excitations 
created by a E are particles, those created by b i _ E are anti-particles, or holes. A 
typical state is 

and such states form a basis for the Fock space of the second quantized theory. 

One can formulate Majorana fermions for a system of this kind if there exists a 
particle-hole symmetry, or, in the relativistic context, a charge conjugation symmetry. 
For example, consider the situation where a matrix T exists such that, for eigenstates 
of H , 

^ E (x) = Vr E (x) (5) 

(This implies that T*r = 1 = rr*.) Then, the particles and holes have identical 
spectra. A Majorana fermion is formed by treating the particle and hole with the 
same energy as a single excitation. The second quantized field operator is 

$(a;,t) = J2 {^E{x)e- iEt l h a E + T^ E {x)e iEt l h a E ^ 

E>0 

This fermion does not have both particles and anti-particles. The ground state |0 > 
is annihilated by a E 

a E \0 >= \/a E 

and a E creates particles, so that the excited states of the system are 

a k a k- a k.l° > 

The field operator is (pseudo-)real in the sense that it obeys 

*(s,t) = r**(a?,t) (6) 
It obeys the anti-commutation relation 

*t(£ t )} =5{x-y) (7) 

To be concrete, in a system of complex fermions where the Hamiltonian such that 
the spectrum has the particle- hole symmetry (J5jl. we could decompose the complex 
fermion into two Majorana fermions by taking the real and imaginary parts, 



$i(s, t) = -= (* (X, t) + T^*(x, t)) 



$ 2 (x,i) = ^(*(x,i)-r**(x,i)) 

Then each of §>\(x,t) and <I>2(x,i) are a Majorana fermion. 

In spite of the beautiful simplicity of this idea, Majorana fermions are not easy to 
come by in nature. One could, for example, decompose the relativistic electron, whose 
wave equation does have a charge-conjugation symmetry, into its real and imaginary 
parts. However, the interaction of the electron with photons is not diagonal in this 
decomposition. The real and imaginary components would be rapidly re-mixed by 
electromagnetic interactions, they cannot be stationary states of the full Hamiltonian 
of quantum electrodynamics. 

One place where we might have better luck is to look for emergent Majorana 
fermions in quantum condensed matter systems. For example, in a superconductor, 
the electromagnetic interactions are effectively screened. Indeed, the Bogoliubov 
quasi-electrons in a superconductor behave like neutral particles. However, even 
there, in an ordinary s-wave superconductor, the anti-particle of a quasi-electron is 
another quasi-electron with opposite spin. Indeed, the quasi-electron operator in an 
s-wave superconductor is the two-component object 

where (|,|) denotes spin up and down. It obeys the charge conjugation condition 

which is not an analog of the Majorana condition in eqn. since it entails both 
conjugation and a flip of the spin. 

In order to find a medium where the quasi-electron is a Majorana fermion, we 
need to consider a superconductor where the condensate has Cooper pairs with the 
same spin, so that the quasi-electron has the form 

f^(x)\ 

Then, quasi-electron is pseudo-real, complex conjugation of its wave- function is equiv- 
alent to multiplying by the matrix T = 

This gives a physical realization of a Majorana fermion. An example of such a super- 
conductor is one with a p-wave condensate, such as Strontium Ruthenate |12j . There, 
the condensate has the form {ip^{x)x x Vip^{x)) (and can in principle have an admix- 
ture of spin down as well). Thus, we see that, in such a material, the quasi-electron 
is a two-component object obeying a Majorana condition. We will make use of this 
example later in this Paper. 




Our particular interest in the following will be in situations where the fermion 
spectrum has mid-gap, or zero energy states. These are well known to lead to in- 
teresting phenomena. Already for complex electrons, mid-gap states give rise to 
fractional quantum numbers |13l I14j . With Majorana fermions, they are known to 
lead to peculiar representations of the anti-commutator algebra which can violate 
basic symmetries |151 116j . Some interesting effects in the context of zero modes on 
cosmic strings have also been examined |17|-|20|. 

To illustrate, let us consider the second quantization of a complex fermion whose 
spectrum has a zero mode, 

H t/j (x) = 
The conjugation symmetry implies that 

If the fermion is complex (not Majorana), the second quantized field has a term 
with the zero mode wave-function and an operator, the first term in the following 
expansion: 



*(x,t) = Vo(x)a + ^E(x)e- iEt ' h a E + ^ E {x)e~ iEt ' h b[ 
Here, a obeys the algebra 



-E 

E>0 E<0 



{a,a f } = l (8) 



and it anti-commutes with all of the other creation and annihilation operators. The 
existence of this zero mode leads to a degeneracy of the fermion spectrum. The 
vacuum state is annihilated by all of the annihilation operators a E and b E . However, 
it must also carry a representation of the algebra © • The minimal representation is 
two-dimensional. There are two vacuum states, (| |>, | |>), which obey 

a E \ T>= = a E \ i> , b E \ |>= = b E \ |> 

and 

a 1 ] i>= | T> , a ] \ |>= 
a\ i>= , a\ |>= | |> 
The entire spectrum has a 2-fold degeneracy, with two towers of excited states, 



and 



°k-°k» 6 k- 6 kl ^> 



having the identical energies ^ • Ei . 

This quantization of the zero mode a was argued by Jackiw and Rebbi |13j to 
lead to states with fractional fermion number. Indeed, the suitably normal ordered 
second quantized number operator 

Q= fdx- ^(x,t),ip(x,t)^ = ^ (a^a E - b^_ E b- E ^j +cJa- \- (9) 

J E>0 



has fractional eigenvalues, for example 



QIT>=+2lt> , Q\i>=-^\i> 

In actuality, the charge operator is defined only up to an overall additive constant. 
However, there does exist a symmetry of the theory, gotten at the second quantized 
level by replacing ^f(x,t) by T^*(x,t). This transformation interchanges particles 
and anti-particles, and is a symmetry of the suitably normal ordered second quantized 
Hamiltonian. It should also flip the sign of Q. It implies that, if there is an eigenstate 
of Q in the system with eigenvalue q, 

Q\q >= q\q > 

then there must exist another eigenstate | — q > in the spectrum of Q with eigenvalue 
-q: 

Q\ - q >= -q\ - q > 

In addition, it is easy to argue that the eigenvalues of Q are space by integers, 
i.e. if qi and q2 are any two eigenvalues, then q\ — q2 =integer. This is essentially 
because the raising and lowering operators for Q are \fd anc i respectively and they 
raise and lower in units of integers. In particular, this implies that 

q — (—q) = 2q = integer 

Thus, the only possibilities are that the entire spectrum of states have integer 
eigenvalues of Q, q =integer, or the entire spectrum of states have half-odd-integer 
eigenvalues q = ^-odd integer. It is easy to see that the operator Q as written in 
© indeed flips sign if we interchange cle <-> &e and a <-» cr and the offset of -1/2 
that appears explicitly there is essential for this transformation to work. This leads 
to the conclusion that, with a single fermion zero mode, the fermion number charge 
is quantized in half-odd-integer units. 

Now, consider what happens for a Majorana fermion with a single zero mode. 3 
In this charge analogous to Q is not defined, so the issue of fractional charge 

is not relevant. But the quantization of the system is still interesting. The second 
quantized operator is 

*(x,t) = M^)a + £ ^ E (x)e- iEt / h a E + £ <M*)e"^/V E 

E>0 E<0 

This fermion contains half of the degrees of freedom of the previous complex one. 
Here, the b E are absent and the zero mode operator is real, a = . 
The creation and annihilation operator algebra is now 

{oe:,c4,} = 8 ee' 

as before, and 

a 2 = 1/2 , {a, a E } = = ja, a^} (10) 
3 We will later construct an explicit example where this precisely this situation occurs. 



A minimal representation can be constructed by denning a vacuum state where 



a E \0 >= for all E > 
Then, we can represent the zero mode by the operator 



a = l(-l)E£>o»W (11) 
V2 



Indeed 
and, since 



a = 



^a E a E \0 >= 

E 

we have 

1 , 

a >= — = > 

The Klein operator, (— l^E>o a E a E ^ anti-commutes with cie and a E . A basis for the 
Hilbert space consists of the vacuum and excited states which are obtained from the 
vacuum by operating creation operators 

a k a k- a kl° > 

These are eigenstates of J2E>o a E a E w hh integer eigenvalues. Thus, in this basis, 
a 2 = 1/2 when operating on each basis vector, and thus the identity operator on the 
whole space. The operator in Ijlljl thus satisfies the algebra Ijlfljl . 

Another, inequivalent representation can be obtained by starting with 

a = — L(_i)E B >o«W (12) 

and a similar construction leads to a Hilbert space whose states are orthogonal the 
one found above. We emphasize here that there are two inequivalent representations 
of the anti-commutator algebra, one where the zero mode operator is represented by 
a in eq. (|TT|) and one where it is represented by a in eq. (|T2|) . Both of these give an 
irreducible representation and the two representations are not related to each other 
by an internal automorphism. 

We observe that these minimal representations of the anti-commutator algebra 
have the property that they break a symmetry of the fermion theory under &(x, t) — > 
— (x, t) , which we shall call "fermion parity" . Fermion parity is a symmetry of the 
linear wave equation even when $(x, t) is a Majorana fermion. At the quantum level, 
fermion parity symmetry leads to a conservation law for the number of fermions 
modulo 2. By this conservation law, any physical process must entail creation or 
destruction of an even number of fermions. For example, if a quantum state is 
initially prepared with an even number of fermions, after any physical process, the 
number should remain even. In operator language, there should exist an operator 
(—1) F which anti-commutes with $(x, t), 

(-l) F $(x,t) + $(f,t)(-l) F = 



and which therefore commutes with the full second quantized Hamiltonian, 



{-\) F H = H{-l) F 

where 

H = J dx^: &(x,t)H $(x,t) : 

However, we see that in the minimal representations of the anti-commutation 
algebras (fTU|) discussed above, in the first representation 

<0|$(x,t)|0>=+-i=V>o(z) 
and in the second representation (|12j) 

<0|$(z,t)|0>=--^o(z) 

In both of these representations, neither the vacuum state, nor any of the excited 
states can be eigenstates of fermion parity, the operator (— 1) . Thus fermion parity 
symmetry is broken by the minimal quantization of this model. 

Fermion parity is a sacred symmetry of physics in four dimensional space-time |21| . 
All fundamental fermions in nature have half-odd-integer spin. A flip in sign of all 
fermion operators can then be realized as a rotation by an angle 2tt. Nature should 
be symmetric under a rotation by 2tt. This means that, if we superpose a state with 
even fermion number and a state with odd fermion number, 

cijeven > +C2|odd > 

no experiment should be devisable, even in principle, to measure the relative sign of 
c\ and C2- In the four dimensional world, unless rotation invariance is broken at a 
the level of fundamental physics, we should always be free to insist that (— 1) F is a 
good symmetry and that we can take all physical states as eigenstates. Of course, 
this applies in four space-time dimensions. The emergent Majorana fermions that we 
want to consider here are embedded in four space-time dimensions. We therefore feel 
free to insist on fermion parity. 

This brings up a contradiction with the previous discussion, where we found that 
fermion parity is necessarily broken by the quantization of the zero mode Majo- 
rana fermion system. The only way to restore the symmetry is to use a reducible 
representation of the anti-commutator algebra. The minimal modification of the 
representation is equivalent to the introduction of another degree of freedom - and 
subsequent use of irreducible representations. The new degree of freedom acts like a 
hidden variable. In the anti-commutator algebra it would be another anti-commuting 
variable /3 which has identical properties to a, 

P 2 = 1/2 

and anti-commutes with all other variables. Then the algebra of a and j3 would have 
a two dimensional representation which we could find by considering the fermionic 



oscillators 



a = -^(a + if3) , a ] = -j= (a - i/3) (13) 
a = -L(a + at) , /3=-^(a-at) (14) 



which obey 



a 2 



, a t2 = , {a, a 1 "} = 1 



We could then find a vacuum state which is annihilated by a, and another state which 
is created from the vacuum by a', 

a\- >= , a)\- >= |+ > 

a\+ >=\- > , at\+ >= 

so that both are eigenstates of (— l)^ and fermion parity is restored. Later we will 
see that the hidden variable (3 can have a physical interpretation. 



2 Degeneracy, tunnelling and teleportation 

In this paper, the most speculative use of Majorana fermions that we shall find is 
for a kind of teleportation by quantum tunnelling. In the context in which quan- 
tum tunnelling is normally studied, a classical object can exist in allowed regions. 
There exist other forbidden regions where it is not allowed to be. Then, quantum 
tunnelling makes use of the fact that, when the particle is quantum mechanical, its 
wave-function does not necessarily go to zero in a classically forbidden region, but 
decays exponentially. That means that it could, in principle, have support on the 
other side of such a region and there is some small probability that an object will be 
found on the other side. This is called tunnelling. 

One might try to make use of quantum tunnelling to transport an object through 
a classically forbidden region. Unfortunately, the exponential decay of the wave- 
function across any classical barrier of appreciable size renders it too small to be of 
any practical use in this regard. A more sophisticated approach would be to create a 
scenario where the wave-function has peaks of appreciable size at spatially separated 
locations, perhaps with a forbidden region in between. This too will fail, but for 
a more sophisticated reason which, since it is related to our later use of Majorana 
fermions, we will outline. Consider, for example the double well potential depicted in 
Fig. ^ If the locations of the minima are well separated and the barrier in between 
them is large, semi-classical reasoning can be applied to this system. Then, the 
ground state of a particle in this potential should indeed have a peak near each of the 
minima, and should be approximately symmetric under interchanging the locations 
of the minima. The typical profile of such a wave-function is drawn in Fig. |2j 

Now, we ask the question. Is this state of use for tunnelling? If this were the 
energy landscape in which a quantum mechanical particle lived, could we, for example, 
populate this ground state by interacting with the system in the vicinity of minimum 1 



/ 2 

Figure 1: A double- well potential. 




and then depopulate the state by interacting with the system near the other minimum 
- 2, effectively teleporting the particle from location 1 to location 2? 

The answer to this question is 'no'. The reason for this answer is degeneracy, 
or approximate degeneracy of the quantum state that we are considering. In such a 
system, when our classical reasoning is good, there must always be a second state, 
perhaps at slightly higher energy but approximately degenerate with the ground state, 
whose wave-function is approximately an anti-symmetric function of the positions of 
the minima. Its typical profile is depicted in Fig. |21 The ground state wave-function 
has the form ifto(x) = ipi(x) + tp2(x) where ipi(x) is localized near minimum 1 and 
ip2(%) is localized near minimum 2. The anti-symmetric state would have the form 

1p a {x) = 1pl{x) - lfj 2 (x)- 

Now, when we interact with the system near minimum 1, while we overlap the 
ground state wave-function, ipo(x), we also overlap ijj a {x) by the same amount. Of 



Figure 3: The almost degenerate state tp a also has two peaks but with differing sig: 



course, the state that we actually populate is a linear combination of the two, 



whose wave-function is entirely localized at the position of the first minimum. The 
particle initially has zero probability of appearing near the second minimum. Our 
attempt at teleportation by tunnelling has been foiled by degeneracy. 

Anytime the Schrodinger equation can be analyzed semi-classically in this way, 
it seems to have a built in protection against the long-ranged behavior that we are 
looking for. 

In this argument, because it is a superposition of two stationary states with slightly 
differing energies, ipi(x) is not a stationary state. It should have a small time de- 
pendence which eventually mixes it with ^(x). But this time dependence mixes it 
slowly, in fact its origin is just the conventional tunnelling amplitude for the particle 
to move from location 1 to location 2 through the barrier in between. 

What we need to find is a quantum system where a quantum state which is well 
isolated from other states in the spectrum can have peaks at different locations. From 
the argument above, it will be difficult to find states of this kind which obey the regu- 
lar Schrodinger equation. Where we will look for such states is in quantum condensed 
matter systems, where electrons, or more properly quasi-electrons, can satisfy equa- 
tions that are very different from the Schrodinger equation. To motivate this, in the 
next Section, we review some of the pictorial arguments for the appearance of frac- 
tionally charged states in polyacetylene. Also, to set the stage for what comes next, 
we discuss what happens if the fermion spectrum of polyacetylene were Majorana, 
rather than complex fermions. 




3 The polyacetylene story 



Before we consider a more quantitative model which will illustrate our point, we pause 
to recall the example of the conducting polymer, polyacetylene. Polyacetylene is a 
hydrocarbon polymer where each Carbon atom bonds with a Hydrogen atom and as 
well forms two strong covalent bonds with neighboring Carbon atoms. The fourth 
valence electron is nominally a conduction electron. However, a Peirls instability 
localizes it into a charge density wave which is effectively a dimer. The result is a 
gap in the electron spectrum at the fermi surface and, without impurities or other 
structures, the material is an insulator. There are two degenerate ground states, 
depending on the direction chosen by the dimerization. We illustrate these as the A 
and B phases in the diagram in Fig. @J In that figure, each line is a covalent bond, 



Phase A 



Phase B 



Figure 4: The two degenerate ground states of polyacetylene. 



using two of the valence electrons of the Carbon atoms. 

The conductivity of doped polyacetylene that is seen by experiments is thought 
to be mostly attributed to the transport of charged solitons along the polyacetylene 
molecules. A soliton in this system is a defect which interpolates between the two 
phases. We have depicted a soliton-anti-soliton pair in Fig. |SJ Note that it can be 



Phase A 



Phase B 



Phase A 



Figure 5: Solitons form phase boundaries. The soliton anti-soliton pair can be created by 
nipping the direction of the bonds between them. 



obtained from one of the ground states by flipping the direction of the bonds that lie 



between the locations of the solitons. Also note that the energy of the system could 
be higher than that of the ground state, since the defects have non-minimal energy 
configurations, but the energy density should be concentrated in the vicinity of the 
solitons. Although it will not be an issue for us, since we are interested in other 
aspects of this system, the solitons turn out to be quite mobile. They also carry 
electric charge, and can thus account for the high conductivity that is attainable in 
polyacetylene. The density of solitons can be controlled by doping. For some original 
literature on polyacetylene, see refs. [22|-[2II1- 

There is a simple argument that shows that a soliton of polyacetylene has half 
of the quantum numbers of an electron In this argument, we will neglect the 

spin of the electron. Thus, for the purpose of our arguments, in figures 0] and each 
bond stands for a single electron, rather than a spin up, spin down pair of electrons. 
Now, consider what happens when we add an electron to phase A, as in Fig. H3 By 




Phase A added electron Phase A 

Figure 6: Phase A with an additional electron. 



flipping the directions of some bonds, we can see that we have created a soliton-anti- 
soliton pair, where each object seems to share half of the added electron. This state 
is depicted in Fig. [Tj 



Phase A 



Phase B 



Phase A 



Figure 7: Beginning with phase A and an additional electron, as shown in Fig. H3 we create 
a soliton-antisoliton pair which seems to share the electron. 



This brings up the question, is the electron really 'split' between the two sites? Or 
does it exist in an entangled state of some sort which has some probability - | - of the 
"whole" electron being located at either site. This question can be made more precise 
by asking about measurement of the electron charge, which is a conserved quantum 
number in this system. If, by further flipping bonds, we separate the solitons to a large 
distance, and then measure the electron charge in the vicinity of one of the solitons, 
is the result of the measurement -e/2? Or does this measurement manage to collapse 
the electron wave-function somehow so that the result is either or -e? In the latter 
case, the average of many measurements might be -e/2, but any single measurement 
would either see a whole electron or no electron at all. The answer to this question 



was found long ago in ref. |3L)l 131 j . The conclusion was that the measurement of the 
electron charge localized near one of the solitons should yield -e/2. Put equivalently, 
the fractional charge of the soliton is a sharp quantum observable. How it manages 
to do this is interesting, and was discussed in ref. We shall review it here. 4 

This issue has recently been reexamined [3B] in conjunction with some ideas about 
entangled electron states in Helium bubbles [15] . 

The electron spectrum in polyacetylene has an electron-hole symmetry. We could 
have created a state with the same energy as the one depicted in Fig. [7|by removing, 
rather than adding an electron, to give a hole which is apparently split between the 
soliton and anti-soliton, as shown in Fig. |HJ 




Phase A 



Phase B 



Phase A 



Figure 8: A soliton- antisoliton pair with a deficit of one electron. 



There are apparently four different states of the soliton-anti-soliton system. There 
are the two overall neutral states, one of which is depicted in Fig. |S] and the other 
obtained by flipping the intermediate bonds in the opposite direction. We could 
also go from one of these states to the other by transporting a whole electron from 
one soliton to the other. The other two states we can obtain by either adding or 
subtracting an electron from one of the ground states and are those that we have 
already discussed in Figs. [7| and |H1 We can identify these charged soliton states in 
the low energy electron spectrum. In the single electron spectrum of polyacetylene 
with a soliton-anti-soliton pair, there are two near-mid-gap states which have small 
positive and negative energies. Thus the low energy electron spectrum has four states, 
a ground state, an electron state, a hole state and an electron-hole state. 

By their quantum numbers, the electron and hole states can be identified with the 
configurations in Figs. [7|and|SJ respectively. The ground state and the electron-hole 
state are neutral and must be formed from linear combinations of the two neutral 
states. Then, in the electron state, the electron wave-function indeed should have 
two peaks, as depicted in Fig. E3 Similarly the hole wave-function also should have 
two peaks, as is depicted in Fig. Detailed analysis shows that one is an even and 
the other is an odd function of relative distance, as shown in the figures. 

The electron wave-function has peaks at two locations. So, we could ask the 
question again: Can we use this system for teleportation? Could we populate the 
electron state by interacting with one of the solitons and subsequently extract the 
electron again, and thereby teleport it, by interacting with the other soliton? To 
understand the answer, which will be 'No!', it is necessary to realize that, once the 



4 For some other literature on this and closely related issues, see refs. |32j-|47|. 
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Figure 9: The electron wave-function. 
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Figure 10: The hole wave-function. 



electric charge is a sharp quantum observable, the electronic states of the solitons are 
disentangled by a local measurement of an observable such as the charge. 

To see how this happens, let us consider a second quantization of this system. 
The electron operator has the form 

ip(x,t) = ip e (x)a + ij)h{x)b* + ... 

where we have identified an electron annihilation operator a for the positive energy 
state and a hole creation operator for the negative energy state. We have neglected 
the time dependence (the energies of the two states are exponentially small in the 
soliton separation). We could as well write 

i/)(x, t) = ifii(x)(a + b ] ) + ip 2 (x)(a - b ] ) + ... (15) 

where ipi(x) = (ip e (x) +iph{x)) has support only in the region of the left-hand 
soliton and ij)2(x) = 4= (fp e ( x ) ~ ^hix)) has support only near the right-hand soliton 
in Figs. Handlini 

Now, if we concentrate on the region near the left-hand soliton, if)(x, t) or ijr(x, t) 
will annihilate or create an excitation using the combination of operators 



' ( a + 6 fN ) , a 1 = — = (a) + b 



o 



s/2\ J ' V2 



Similarly, if we concentrate on the region around the right-hand soliton, excitations 
are created and annihilated using 



The set of operators (a, a', [3, /?') are a Bogoliubov transformation of the creation 
and annihilation operators (a, ,b,b^). This transformation does not violate fermion 
number - it superposes operators with the same fermion number. Further, the ex- 
citations that the new operators create or annihilate are entirely localized on one or 
the other of the solitons. 

Thus, again, we do not have a process whereby an electron or hole state which 
has two peaks can be populated by interacting with the system in the vicinity of one 
of the peaks. We have failed to find teleportation. Instead we have found fractional 
charge. The charge density integrated over the vicinity of one of the solitons turns 
out to be 

Q = — e (afia — l/2j + charge of electrons — charge of holes 

which indeed has half-odd-integer eigenvalues. This Bogoliubov transformation, as a 
mechanism for disentangling the charge quantum numbers of the solitons was origi- 
nally found in ref. |31j . 

In eq. (|15|). we ignored the small time dependence of the near mid- gap states. At 
this point, the reader might wonder if the disentanglement of the soliton and anti- 
soliton charges that we find by the Bogoliubov transformation would not be undone 
by this time variation. Indeed, it would be, eventually. However the time scale is 
given by the inverse of the energy gap and is therefore exponentially large in the 
distance L between the soliton and anti-soliton, T ~ m~ 1 e mL , where m is the energy 
gap. This is roughly the time for quantum mechanical tunnelling between the solitons 
assuming an energy barrier of height the energy gap extending over distance L. For 
macroscopic L this time T should be very large. 

What has prevented teleportation in this second example is again a degeneracy, 
this time a slightly more subtle one since, even though the electron and hole state 
have identical energies, they have opposite signs of charge. Avoiding teleportation has 
led to fractional charge. It has done this by a hybridization, at the second quantized 
level, of the propensity of the electron field operator to create an electron and to 
annihilate a hole in a local state. 

Now, imagine that, rather than complex electrons, polyacetylene had Majorana 
fermions which would be obtained by identifying the particle and hole states as the 
same excitations. (Here, we are ignoring the obvious disaster that this scenario would 
lead to in chemistry.) Then, in eqn. we would have to identify a = b and 

^MajOM) = ipi(x)(a + a" 1 ") + ip 2 {x){a - a)) + ... (16) 

Now, a + a"!" cannot be an annihilation operator, in fact 



{a + a)f = 1 



It is similar to the single zero mode operator "a" that we found for a Majorana 
fermion in the Eq. (|14|) . In fact, the other combination — a ) now plays the role 

of u j3" , the "hidden variable" . Its purpose in our previous discussion was to provide 
a quantization which did not violate fermion parity. Here, this hidden variable is 
just the fermion zero mode living on the far-away anti-soliton. We could choose the 
fermion parity conserving quantization by using the states (|— >, |+ >) defined by 

a\- >= , a)\- >= |+ > 

a\+ >= |- > , a)\+ >= 

which can be eigenstates of (—1) F . In these states, the expectation value of the 
fermion operator vanishes, for example < 0\ipMa](x, t)\0 >= 0. However, the two 
solitons are invariably entangled. There is now no conserved fermion number that we 
can use to measure this entanglement, but there are other effects which we will discuss 
in later sections once we have made the present reasoning more solid by discussing 
it in the context of a field theoretical model and them formulated a more realistic 
model with emergent Majorana fermions. 

4 Relativistic Majorana fermions in a soliton 
background 

Single-particle states that are in some sense isolated are well known to occur for 
Dirac equations, particularly when interacting with various topologically non-trivial 
background fields such as solitons, monopoles and instantons. The consequences of 
fermion zero modes such as chiral anomalies [SU] and fractional fermion number |13| . 
|14j are well known. 

The polyacetylene example, in the context of discussions of fractional charge, 
that we used in the previous Section is a well-known example of this. In polyacety- 
lene, the low energy electron spectrum can be approximately described by the Dirac 
equation |221 151 j and the solitons which we discussed using pictures have a mathe- 
matical description as soliton-like configurations of a scalar field which couples to the 
Dirac equation. In this Section, we will make the analysis of the previous Section 
more quantitative by considering the problem of a 1+1-dimensional relativistic Dirac 
equation coupled to a soliton background field and a soliton-anti-soliton pair. 

Consider, for example, the simple one-dimensional model with Dirac equation 

[%>fd li + <j>(x)]ij>(x,t) = (17) 
The Dirac gamma-matrices obey the algebra 

{7",7"} = 2<r 

where g^ v = is the (inverse of the) metric of two dimensional space-time. 

This describes a fermion moving in one dimension and interacting with a scalar 
field <fi{x) which we shall take to have a soliton-anti-soliton profile. For the purposes 



of this discussion, we take the ideal case of a step-function soliton located at position 
x = and a step- function anti-soliton located at x = L, 



60 x < , x > L 
00 < x < L 



(18) 



We will assume that the solitons are very massive, so they do not recoil when, for 
example, fermions scatter from them. 
If we take 

1>{x,t) =^ E (x)e- lEt 



and choose an appropriate basis for the Dirac gamma-matrices, the Dirac equation 
becomes 

- 4>{x) J \v E (x), 
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u E {x) 
v E (x) 



(19) 



This equation has a particle-hole symmetry 

ip- E {x) = i)* E (x) 

It is easy to show that it has exactly two bound states. One is a state with small 
positive energy and the other is the associated hole state with a small negative energy. 
The wave-functions 



(20) 
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Q e~^ x + e*°i x - L ) + 0{er^ L ) < x < L (21) 
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(22) 



e -<t> x + Q e ^{x-L) + 0( e -^ L ) < x < L (23) 
^ e ML-x) + o( e -*oi) L < X 



where, sufficient for our purposes, we give only the large L asymptotics - corrections 
to all quantities are of order e~^° L . Note that ip-(x) is indeed related to ip+(x) by 
ip-(x) = ip+(x). 

These states have energy well separated from the rest of the spectrum, which is 
continuous and begins at E = ±<f>o. The energies are also exponentially close to zero 
as the separation L is large. Furthermore, each wave-function has two peaks, one near 
x = and one near x = L. They have identical profile near x = and they differ by 



a minus sign near x = L. This is the same feature of the electron and hole states that 
we claimed for the polyacetylene soliton-anti-soliton system in the previous Section. 
The second quantized Dirac field now has the form 



il>(x, t) = ^ + {x)e- iE+t a + ^* + {x)e iE+t tf + ... 



(24) 



When L is large, one can consider a second set of almost stationary states which 
are the superpositions 
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< x < L (26) 
L < x 



which has most of its support near x = and 
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which has most of its support near x = L. 

In terms of these wave-functions, which are localized at the sites of the solitons, 



Y>(x, t) = Y> (x, t)-j= (a + 6 f ) + Vl(^, *)^=t (-a + fot ) + • • • 
We could now consider the creation and annihilation operators 



(29) 



a = i=(a + 6t) , a t = J_( a t +6 ) 



1 / 
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By interacting with the system at x = 0, we could as well be dropping the fermion 
into the state ipo, which is localized there and which has exponentially vanishing 
probability of occurring at x = L (until sin Eot becomes appreciable, which is just 
the usual estimate of tunnelling time through a barrier of height (fro and width L). 

It might seem bizarre that, if we begin with the system in its ground state when 
L is small, then adiabatically increase L that we would not simply end up with the 
original ground state that has ip-(x) populated, ip+(x) empty In fact, this is a 



possibility. However, as we have argued in the polyacetylene example in the previous 
Section, as L — » oo, the result is an entangled state of (appropriately defined |3U1I3T] 1 
fermion number. If we begin with the original ground state, measurement of the 
fermion number which is localized in the vicinity of one of the solitons will collapse 
the wave-function to one where the fermion, rather than occupying the negative 
energy state tp-, occupies either the state ipo or the state ipL which are localized at 
x = or x = L, respectively. As seen from the vicinity of each soliton, these are 
identical to the Jackiw-Rebbi states |13j of the fermion in a single soliton background, 
which have fermion number These states are time-dependent, but again, just 
as in the polyacetylene example, the time scale for charge fluctuations is just the 
tunnelling time for a particle to go between the locations of the solitons. 

What about teleportation? Now, our dumping a fermion into the bound state, if 
performed near x = would populate the state tpo(x), rather than ip + (x), as all local 
operators would couple only to this state. It would have appreciable probability of 
appearing at x = L only after a time over order Eq 1 ~ <pQ l e^° L . 

The situation is somewhat different if we assume that the fermion is a Majorana 
fermion. The Hamiltonian of a Majorana fermion must have a symmetry which 
maps positive energy states onto negative energy states. In the case of (|19|) . we have 
ip-E{x) = i(j* e {x). Then, a fermion and an anti- fermion have the same spectrum, and 
we can identify them as the same particle. 

Now, for the Majorana fermion, the pair of wave- functions ip+(x) and ip-{x) 
correspond to the same quantum state which can be either occupied or empty. (We 
can arbitrarily assign fermion parities (— 1) F = — 1 for the unoccupied state and 
(—1) F = 1 for the occupied state, although +i and —i might be more symmetric). In 
this case, the states ipo and ipL are wave-functions for superpositions of the occupied 
and unoccupied states - they do not have definite fermion parity. 

If we begin with the system where the quantum state is an eigenstate of fermion 
parity and we by some process dump a fermion into the bound state near x = 0, its 
wave-function automatically has a second peak at x = L and it could in principle be 
extracted there. This defines what we mean by "teleportation" . 

If we concentrate on the region near x = and we are unaware of the region near 
x = L, depending on the quantization, this teleportation will appear as either viola- 
tion of conservation of fermion number mod 2 or the existence of a hidden variable 
in the local theory. 

5 P-wave superconductor model and Andreev 
states 

Of course, the fermions in polyacetylene are not Majorana, they are electrons with 
complex wave- functions. The place to look for emergent Majorana fermions in nature 
is in superconductivity. Here we shall formulate a model whose basic excitations are 
Majorana fermions. We will do this by using contact with a p-wave superconductor to 
violate the conservation of total charge, leaving behind conservation of charge modulo 
2. In such an environment, the real and imaginary parts of the electron can have 



different dynamics and the electron is essentially split into two Majorana fermions. 
They can further be coupled to soliton-like objects, in this case the boundaries of the 
space, in such a way that only one of the Majorana fermions has zero modes. Then, 
the scenario that we have been looking for, an isolated single-particle state, can be 
found. 

In these materials, mid-gap bound states, called Andreev states, are a common 
occurrence. They typically live at surface of the superconductor [^2j. In our case, 
these will be Majorana zero modes. 

Majorana zero modes of the type that we are discussing are also known to be 
bound to vortices in p-wave superconductors where they have the remarkable effect 
of giving vortices non-Abelian fractional statistics [53], [51]. For concreteness we will 
consider a slightly simpler model one-dimensional model that was originally discussed 
by Kitaev [2] in the context of fermionic quantum computation. 

We shall consider a quantum wire embedded in a bulk P-wave superconductor as 
is depicted in Fig. ^3 




Figure 11: A quantum wire embedded in a bulk P-wave superconductor. 



We shall assume that the wire has a single channel. We shall also assume that 
the dynamics of electrons in the wire are adequately described by a one-dimensional 
tight-binding model. We will ignore the spin degree of freedom of the electron. The 
phenomenon that we will find is to a first approximation spin-independent. 

We will assume that the coupling to the neighboring p-wave superconductor is 
weak and its net effect is to give electrons the possibility of entering and leaving the 
wire in pairs by creating or destroying a p-wave cooper pair in the bulk. To describe 
the electrons, we will use the Hamiltonian 

( t t* A A* \ 

H = ^2 ( -a ] n+l a n + — a j n a n+ i + —a ] n+l a) n + —a n a n+ i + ^a) n a n J (30) 

n=l V - - - / 

Sites on the quantum wire are labelled by n = 1,2, ...,L. The operators a n and a} n 
annihilate and create an electron at site n. They obey the anti-commutator algebra 

I "' !'} = <W (31) 

The first terms in the Hamiltonian, with coefficients t and t* are the contribution 
to the energy of the hopping of electrons between neighboring sites. The second pair of 
terms, with A and A*, arise from the presence of the super-conducting environment. 
They describe the amplitude for a pair of electrons to leave or enter the wire from 
the environment. It is assumed that they can do this as a Cooper pair when they 
are located on neighboring sites. This is effectively an assumption about the size and 
coherence of the cooper pairs in the superconductor. Even if it were not accurate, 



the smaller next-to-nearest neighbor, etc. terms that would arise could be taken into 
account and would not change our result significantly. The last term is the chemical 
potential, the energy of an electron sitting on a site of the wire. We shall assume 
a reasonable hierarchy of the parameters, that the amplitude for hopping along the 
wire is somewhat larger than hopping to and from the bulk, \t\ > |A|, and that 
the chemical potential is close enough to zero that the electron band has substantial 
filling, \fi\ < \t\. 



5.1 Spectrum of single-particle states 

Let us discuss the spectrum of the single-particle states in the many-body theory 
described by the Hamiltonian (J3UJ). If t = |i|e 4<?i and A = |A|e 2 * e , by redefining 
a n — > e % ^ +e ^a n for n odd and a n — > e l ^~ 9 ^a n for n even, we remove the complex 
phases of t and A, which we can henceforth assume to be positive real numbers. The 
equation of motion for the fermion wave-function is gotten by taking the commutator 
of its operator with the Hamiltonian (J3(J|). 

iha n = [a n ,H] 

for which we get 

d t A / + +\ 

~dt an = 2 ^ n+1 + an ~^ ~ 1 \ n+1 ~ n ~ x J + ' ' 

for the sites n = 2, L — 1. 

Because we are using open boundary conditions - the chain simply ends at n = 1 
and n = L, the equations for and ^a-L differ from (|32j) by missing terms. When 
we solve ()32j) as a wave equation, it will be convenient to deal with this by extending 
the chain by one site in each direction and then eliminating the extra sites by imposing 
Dirichlet boundary conditions, 

a (t) = , a L+1 (t) = 

With these conditions, (|32|) describes the dynamics for all n = 1,2, L. 

Now, it is most efficient to decompose the electron into real and imaginary parts, 
a n = b n + ic n , and assemble them into a spinor 

*• " (t) (33) 

Note that this spinor obeys the Majorana condition 

1> n = p n (34) 
The equation for the wave-function is 

. H £ t ^J^ + V i(t + A )J^ +1 + V o i^-A);^- 1 - 

(35) 



In order to solve the equation, we will make the ansatz 



The Majorana condition for energy eigenstates is 

CM = ^n(-w) 
We will normalize the wave-functions with the condition 



(36) 
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Since the equation and boundary conditions are linear, we can further make the 
ansatz that the wave-functions are superpositions of plane waves, 



n ( U>{() 

v(C) 



(37) 



Then, the difference equation (|35jl becomes 
The equation for the wave-function is 



£t(C + i/0-3 A K-i/0 + M 

it(C + i/c)-iA(c-i/C) + /i. 
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which has a solution when the frequencies obey the dispersion relation 



|t(C + l/C)+A* 



l 2 



^A(c-i/C) 



-i 2 



(38) 
(39) 



For a given real value of oj, there are generally four wave- vectors which satisfy this 
dispersion relation, 

To find a solution of the wave equation, we must take superpositions of the four 
solutions of (|3*5|) with each of these four wave-vectors. Then we must adjust the 
four coefficients of the superposition in order to satisfy the four boundary conditions. 
(Remember that the boundary conditions are for spinors, so there are four boundary 
conditions in total.) Three of the boundary conditions can be solved by adjusting the 
coefficients in the superposition. The fourth superposition coefficient can eventually 
be determined up to phases by normalizing the wave-function. The fourth boundary 
condition, which has yet to be satisfied, then gives a condition that the wave-vector 
must obey. Plugging the resulting wave-vector back into the dispersion relation (|3*9"|) 
then gives the allowed energy eigenvalue. This gives an algorithm for finding the 
energies, the allowed wave- vectors (which are i In £ and are generally complex) and 
the wave- functions. 

When L is large, the solutions are of two kinds. One are to a good approximation 
continuum states, where £ = e lk and the continuum spectrum is 



to(k) = ±J [t cos k + fj] + A 2 sin 2 k 



with k £ (— 7r,7r] (it is quantized approximately as k = 2tt ■ integer/(L + 1) which 
becomes a continuum when L — * oo). This spectrum has an energy gap. The point 
of closest approach of the positive and negative energy bands occurs when cos k = 

—tfi/(t 2 — A 2 ) and the gap is E gap = 2Ay / ^-^^^-. We will assume that this gap 
is significant, so that the mid-gap states that we will discuss next are indeed well 
isolated. 

The other modes in the spectrum are a pair of mid-gap states. When L is large, 
these states have energies that are exponentially small in L, one is positive, one is 
negative and they have equal magnitudes. In the following, we will solve for the 
spectrum of these mid-gap states in the approximation where effects that are expo- 
nentially small in L are neglected. 

We begin with an un-normalized spinor 
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w c * + i/ c * ) + !a(c*-i/c*) + // 



We will solve the boundary condition for the mid-gap state in the limit where L is 
large. There, we expect the solution to be very close to u = 0, for which we then 
need a wave- vector which solves i(£ + 1/C) + 2/U = — A(£ — 1/C)- Then, to a first 
approximation, the terms with A and C are absent from (j4()|) and we must choose 
the B = — (( — 1/C)/(C* — 1/C*) i n order to satisfy the boundary condition at n = 0. 
Since 
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/i 



2(t + A) 



so that (X* = |jp2f < 1, this gives a wave- function which is maximal at n = 1 and 
which decays exponentially as n increases. This would indeed be the solution for the 
mid-gap state on the half-line when L — > oo. When L is finite, rather than infinite, in 
order to satisfy the boundary condition at n = L + l we must include an amplitude for 
the growing solution. It can be obtained from the decaying one by simply replacing 
n by L + 1 — n and multiplying the spinor by a 2 . Thus, to a good approximation the 
mid-gap solution is 
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This wave-function has infinitesimal positive energy. The wave-function with in- 



finitesimal negative energy is given by 
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We will abbreviate these by naming the function 

I a t 2 -^ {-^+ lv w^^y-{-^-zyt 2 -A 2 -^y 

0n " 'V Ytt 2 -A 2 -fi 2 (t + A)™ 

(41) 

where = and we have normalized to 



(42) 



The function 4> n has maximum magnitude at n = 1 and it decays exponentially as n 
increases. We shall use the notation 

V>+ = K (j) - 0L + l-n (o) (43) 



C=^(j)+0L + l-nQ (44) 

We have normalized the spinors so that 

n 

Note that, these wave-functions satisfy the Majorana condition ip~ = As 
expected, they have support near n = 1 and n = L and are exponentially small in 
the interior of the quantum wire, far from the boundaries. Further, we have adjusted 
phases so that the wave- functions are identical in profile in the region near n = 1. 
Then, we expect that they differ in sign in the region near n = L and we confirm 
from that above that this is so. Also, note that they are complex. To form the real, 
Majorana spinor, we must superpose them with a creation and annihilation operator, 

4>n(t) = ipn e ~ luJt(1 + V'n e lw *a^ + non — zero energy states (46) 

Here a and a) are the annihilation and creation operators for the mid-gap state and 
uj is their exponentially small energy. Ignoring the energy, we can also write this 
operator as 

i>n{t) = 4>n (^j (a + a f ) + <f>L+l-n (^j \{d ~ <J) + . . . 



The first term on the right-hand side has support near n = and decays exponentially 
as n increases from 1. The second term has support near n = L and decays expo- 
nentially as n decreases from L. They each multiply the operators a = A= ((a + a*} 

and (3 = (a — a'), respectively. These are analogous to the operators which we 

introduced on Section 1. a and a) must have the anti-commutator 

{a,a f } = 1 

which has a two-dimensional representation, the states |— > and |+ > of Section 1 
which we copy here for the reader's convenience, 

a\- >= , a)\- >= |+ > 

a\+ >= |- > , <J\+ >= 

All other excited states of the system are created by operating creation operators for 
the other, non-zero energy excitations. Remember that it is the states |+ > and |— > 
which we expect to be eigenstates of fermion parity, (—1)^. 



5.2 Second quantized electron operator 

Now, we recall that the upper and lower components of the spinor ip n (t) that we 
discussed in the previous subsection are simply the real and imaginary components 
of the electron field operator, which we can now reconstruct, 

a n (t) = ^L+i-rX (a ~ a f ) + i<Pn (a + a)^j + ... (47) 

This is now a complex operator, but its real and imaginary parts have support at 
opposite ends of the quantum wire. The part of the operator which has not been 
written, and is indicated by dots in (|47|). are superpositions of creation and annihi- 
lation operators for continuum states. All such states have energies above the gap 
and extended, plane-wave-like wave- functions. Note that now that the phase symme- 
try of the system has been broken by coupling to the superconductor, the real and 
imaginary parts of the electron operator will generally have different properties. This 
interesting fact will not concern us in the following and we will focus on the mid-gap, 
or zero mode part of the electron operator. 

Note, now, if we operate with any local operator in the vicinity of n = 1, the 
electon operator acts as if it were composed of the combination of zero mode operators 
(a + a'). As we have discussed before, this operator squares to a constant. There 
cannot be any states that it annihilates. Thus, operating it on any state of the system, 
in the region where the zero mode wavefunction has support, will have an effect. What 
it does is flip the state from |— > to |+ >. Since it is a hermitian operator, it is 
possible to diagonalize it, the states ^(l — > +1+ >) an d > +1+ >) are i* s 

eigenvectors. However, these eigenvectors are not eigenstates of fermion parity. 



6 Long ranged correlations of electrons 

What about teleportation? Let us imagine that we begin with the system in one of 
its ground states, say | — > and inject an electron so that at time T = it is resting 
at site #1. This means, we being with the state a\\— >, where, as we recall, a\ is 
the creation operator for an electron at site #1. 

We then ask what is the quantum transition amplitude for the transition, after 
a time T has elapsed, of this state to one with the electron located at position #L. 
The final quantum state is a} L \ — >. The amplitude is given by 

A =< -| a L e iHT a\ \- >= |0?| 2 + (T and L-dependent) (48) 

The T- and L-dependent parts of this matrix element represent the usual propagation 
via excited quasi-electrons which must travel across the wire. The first term is non- 
zero and is T and L-independent. By 'teleportation', we are referring to this part of 
the amplitude. Here, we can evaluate the amplitude explicitly. It is 

/2A\ ft 2 - A 2 -^\ 
A ™ = {-) { ( t + A)* J 

Which can be appreciable, in the 10-30 percent range, for a surprisingly wide choice 
of parameters. 

However, the teleportation probability is the square of this amplitude, which 
is somewhat smaller. We could ask a more sophisticated question: What is the 
probability that the electron, once injected at n = 1 appears anywhere within the 
exponential range of the zero mode wave-function at n = L. This probability would 
be given by 

PTcl = J>"| 2 ^ 

n 

This is what we shall call the "teleportation probability". Again, for a range of 
parameters t, A and fx, it can be appreciable. 

7 Discussion 

The apparently instantaneous propagation of an electron would seem to be a potential 
violation of Einstein causality, since in principle a message could be sent at a speed 
faster than that of light. 

Let us review the nature of the system that we have constructed. Once the 
quantum wire - p-wave superconductor system is prepared, the extended Majorana 
state of the electron is already there, ready for use. The system has a two-fold 
degeneracy: at low energy, there are two states |— > and |+ >. These are not normal 
quantum states in that they differ by a quantum number which we would like to 
preserve, fermion parity (—1)^. 

Thus, if we do not allow superpositions of these states, this is effectively a classical 
bit, like a classical switch that can either be OFF or ON, the wave-function can be 
in one state or the other. 



(49) 



The system moves from OFF to ON by absorbing or emitting an electron in a way 
that flips the vacuum from one state to the other. This should occur somewhere in the 
vicinity of the ends of the wire, where the zero mode wave-functions have support. It 
can move back from ON to OFF by the identical process, again absorbing or emitting 
a single electron. 

This leads to the rather drastic conclusion that there could be super-luminal 
transfer of information in this system. One would need only to prepare the system in 
one of its ground states, with a sender sitting at 1 and a receiver sitting at L. Either 
ground state is sufficient and neither the sender nor the observer needs to know which 
it is. All the receiver has to do is wait for an electron to arrive. If it arrives with 
energy at or above the electron energy gap, he or she can conclude that it propagated 
normally and was sent at some time in the past. However, if it arrives at very low 
energy, he or she knows that it tunnelled and that it was sent by the sender at that 
instant. This is seemingly an instant transfer of information over a finite distance. 

There is a obvious way out of this, but it means giving up the fermion parity 
symmetry that has until now been sacred. If we allow superpositions of the states 
|— > and |+ > which have even and odd fermion number, then the degenerate ground 
states are a quantum rather than classical two-level system, there are two states and 
any superpositions are allowed. Now, in this system, it is easy to prepare states where 
an electron can spontaneously appear or disappear. Take, for example an eigenstate 
of the operators that we called a and (5. In their eigenstates, ^= (|— > ±| + >), the 
electron operator has an expectation value < a n (t) >= ±(j) n ± <pL+i-n- It would thus 
have an amplitude for simply vanishing or appearing spontaneously. 

Then, when the observer at L detects the arrival of a low energy electron he or 
she cannot distinguish one which was sent from the other side of the wire from one 
which is spontaneously created. This restores Einstein causality at the expense of our 
having to admit states onto physics which are not eigenstates of fermion parity. There 
is the further question of whether such states are consistent with three dimensional 
physics. 

Fermion number mod 2 is an important conservation law in three dimensional 
physics |55j . Even though the quantum wire that we have discussed is one-dimensional, 
it is embedded in three dimensional space and the electrons that we are discussing 
are spin-^ particles in three dimensional space. This means that their wave-functions 
individually change sign under a rotation by 2ir. More importantly, a state with odd 
fermion number must change sign under a rotation by 2tt whereas a state with even 
number should remain unchanged. A rotation by angle 2tt should not affect physics. 
Thus, the relative sign of even and odd fermion number states should not have any 
physical consequences. 

If we did allow a superposition of the two states, they would form a single qubit. 
We could parameterize the state-vector by a point on the Bloch sphere (9, <f>) where 
the state is 

9 9 
\0,<p> = cos-|- > +e itp sin-\+ > (51) 

Points on the two-dimensional unit sphere are specified by the unit vectors 



h = (sin 9 cos </?, sin 9 sin <p, cos 9) 



and < 9 < tt, —tt < ip < tt. However, as we have argued, the relative sign of the 
two states should not be an observable. Then the set of "physical states" of the qubit 
would be the Bloch sphere with a further identification 

if ~ if + TT (52) 

Of course, this identification is allowed only if there are no experiments, even in 
principle, which could measure the relative sign of the two states in the superposition. 
Normally, one could measure that sign by an interference experiment. 

For example, we could attempt to observe the relative sign by examining interfer- 
ence between the electron which arrives by tunnelling and the one which arrives by 
conventional transport. However, the teleportation amplitude in the state eft 

< 6,ip\ cll e rHT a\ \9,ip >= cos 6 [teleportation] + [transport] (53) 

The teleportation amplitude is diminished by a factor of cos 9 whereas the transport 
amplitude is unchanged. One can make the teleportation amplitude vanish by ad- 
justing 9 = tt/2. However, the relative amplitude cannot be used to measure the 
relative sign of the two components of the wave-function. 
There is an amplitude for an electron to vanish, 

< 9, p\ e lHT a\ \9,(p >~ i sin 9 cos ip ■ 4>i (54) 

and to appear spontaneously 

<9,tp\aLe lHT \0,ip >~ sin#sin</? • a (55) 

As we expect, the latter two amplitudes change sign when we put ip — > (p + ir. Actual 
quantum observables are probabilities which are the modulus squares of amplitudes. 
They are also insensitive to the relative sign of the two parts of the wave-function. 

The above probability amplitudes do not offer a way to distinguish the quantum 
states with ip and <p + 7r. At this point, we have not ruled out, but also we have not 
devised an experiment by which they could be distinguished. Indeed, if there is no 
such experiment, we are free to cut the Bloch sphere in half by the identification (|52ft 
and the ground states would form this peculiar qubit. Teleportation still happens, 
but so does the spontaneous disappearance or appearance of a single electron and the 
contradiction with Einstein locality is removed. 

We cannot exclude the possibility that the effect that we have been discussing 
could be interfered with by the superconductor which the quantum wire is in contact 
with. Here, we have assumed that it acts as a simple bath which supplies and absorbs 
Cooper pairs but is otherwise innocuous. We cannot rule out that it also has exotic 
states that should be included in the picture. 
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